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Abstract 

In this paper, we develop a viscosity method for Homogenization of Nonhnear 
Parabohc Equations constrained by highly oscillating obstacles or Dirichlet 
data in perforated domains. The Dirichlet data on the perforated domain can 
be considered as a constraint or an obstacle. Homogenization of nonlinear 
eigen value problems has been also considered to control the degeneracy of 
the Porous medium equation in perforated domains. For the simplicity, we 
consider obstacles that consist of cylindrical columns distributed periodically 
and perforated domains with punctured balls. If the decay rate of the capac- 
ity of columns or the capacity of punctured ball is too high or too small, the 
limit of Me will converge to trivial solutions. The critical decay rates of having 
nontrivial solution are obtained with the construction of barriers. We also 
show the limit of satisfies a homogenized equation with a term showing 
the effect of the highly oscillating obstacles or perforated domain in viscosity 
sense. 
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1. introduction 

This paper concerns on the homogenization of nonlinear parabolic equa- 
tions in perforated domains. Many physical models arising in the media with 
a periodic structure will have solutions with oscillations in the micro scale. 
The periodicity of the oscillation denoted by e is much smaller compared to 
the size of the sample in the media having macro scale. The presence of slow 
and fast varying variables in the solution is the main obstacle on the way 
of numerical investigation in periodic media. It is reasonable to find asymp- 
totic analysis of solutions as e goes to zero and to study the macroscopic or 
averaged description. In the mathematical point of view, the partial differ- 
ential equation denoted by may have oscillation coefficients and even the 
domains, Q^, will have periodic structure like a perforated domain. So for 
each of e > 0,we have solutions satisfying 



I/fMf = in Q 



with an appropriate boundary condition. It is an important step to find the 
sense of convergence of to a limit u and the equation called Homogenized 
Equation 

Lu = mVt 
satisfied by u. Such process is called Homogenization. 



Large number of literatures on this topic can be found in BLPj], JKO |. 



And various notion of con verge nces have been introduced, for example P- 



convergence of DeGiorgi, [DG|, G-convergence of Spagnolo ,[S|, and H- 
convergence of Tartar, Tl|. Two-scale asymptotic expansion method has 
been used to find L formally and justified by the energy method of Tartar. 
He was able to pass the limit through compensated compactness due to a 
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particular choice of oscillating test function For the periodic structure, 
two-scale convergence was introduced by Nguetseng, and Allaire, 0, 
which provides the convergence of u^{x) to a two-scale limit uo{x,y) in self- 



contained fashion. And recent 



y viscosity meth od for homogenization has 

Nonv ariational prob- 



Ca 



CLl| 



CL2|. They observe 



been developed by Evans, Evl[ | and Caffarelli, 
lems in Homogenization has been considered in 
that the homogenization of some parabolic flows could be very different from 
the homogenization process by Energy method. For example, there could be 
multiple solutions in reaction diffusion equations. It is noticeable that the 
parabolic flows with initial data larger than largest viscosity elliptic solution 
will never cross the stationary viscosity solution and that the homogeniza- 
tion will happen away from a station ary solution achieved by minimizing the 
corresponding energy, CLM! | CLM1| . And the viscosity method has been ap- 
plied to the h omogenization of nonlinear partial differential equations with 
random data 



CSWj . |CLM2| . 



Now let us introduce an example of parabolic equations in perforated 
domains. Set f2 be a bounded connected subset of M"" with smooth boundary. 
We are going to obtain a perforated domain. For each e > 0, we cover M" by 
cubes Um<£tZ"Qm "where a cube Ql^ is centered at m and is of the size e. Then 
from each cube, Ql^, we remove a ball, Ta^m), of radius having the same 
center of the cube Q^. Then we can produce a domain that is perforated by 
spherical identical holes. Let 

Rl := R^\Ta, 

and 

Qa, ■.= nnRi = n\Ta,, 
QT,a, ■.= ^a,y< (o,T]. 

Now we are going to construct the highly oscillating obstacles. Let us consider 
a smooth function (f{x, t) in Qt = x (0, T] which is negative on the lateral 
boundary OiQt, i.e. </) < on diQr and positive in some region of Qt- Highly 
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oscillating obstacle ip^{x,t) is zero in fla^ and f{x,t) on each hole 

if gT,^ X (0,T] 

Otherwise. 

Then (f^{x,t) will oscillate more rapidly between and (p{x,t) as e goes to 
zero. 

We can consider the standard obstacle problem asking the least viscosity 
super-solution of Heat operator above the given oscillating obstacle: find the 
smallest viscosity super-solution u^{x,t) such that 



H[u] = An, -ut<0 inQr {=nx (0, T]) 

u,{x,t) = ondiQr (=5fix(0,r]) 
Ue{x,t) > (pe{x,t) in Qt 

u^{x, 0) = g{x) onQ X {0} 

where g{x) > ip{x,0), ipe{x,t) < on diQx and ip^ is positive in some region 
of Q t- The concept of viscosity solution and its regularity can be found at 



CC|. 



We are interested in the limit of the as e goes to zero. Then there are 
three possible cases. First, if the decay rate a, of the radius of column is too 
high w.r.t. e, the limit solution will not notice the existence of the obstacle. 
Hence it will satisfy the Heat equation without any obstacle. Second, on the 
contrary, if the decay rate a, is too slow, the limit solution will be influenced 
fully by the existence of the obstacle and then become a solution of the 
obstacle problem with the obstacle (p{x). We are interested in the third case 
when the decay rate a, is critical so that the limit solution will have partial 
influence from the obstacle. Then we are able to show that there is a limiting 
configuration that becomes a solution for an operator which has the original 
operator, i.e. Heat operator, and an additional term that comes from the 
influence of the oscillating obstacles. Naturally we ask what is the critical 
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rate a* of the size of the obstacle so that there is non-trivial limit u{x) of 
u^{x) in the last case and what is the homogenized equation satisfied by the 
limit function u. 

The elliptic variational inequalities wi th h ighly oscillating obstacles were 
first studied by Carb one ari d Colombini, and developed by DeGiorgi, 

Dal Maso and Longo, |dDi| . Dal Maso and Longo, |dl|, Dal Maso, |di| pj. 



H. Attouch and C. Picard, 



APj ] , in more ge nera. 



was considered b y Cio ranescu and Murat, CM| |. The other useful references 



context. The energy method 



can be found in CM| . The method of scale-convergence was adopted by J. 
Casado-Di'az for nonlinear equation of p-Laplacian type in perforated domain, 
and the parabolic version was studied by A.K. Nandakumaran and M. Rajesh 



NR| . They considered the degeneracy that is closed to parabolic p-Laplacian 



type and that doesn't include the Porous Medium Equation type. L. Baffico, 
C. Conca, and M. Rajesh considered homogenization of eigen value problems 
in perforated domain for the nonlinear equation of p-Laplacian type, BCR |. 

The obstacle problems for linear or nonlinear equation of the divergence 
;ype has been studied by many authors and the reference can be founded in 
Frl |. The viscosity method for the obstacle problem of nonlinear equation of 
non- divergence type was s tudied by the author jLl| . L2|. 



Caffarelli and Lee 



ClJ] develop a viscosity method for the obstacle prob- 



lem for Harmonic operator with highly oscillating obstacles. This viscosity 
method is also improved into a fully nonlinear uniformly elliptic operator 
homogeneous of degree one. 

The homogenization of highly oscillating obstacles for the Heat equation 
has been extended to the fully nonlinear equatio ns o f non- divergence type. 
This part is a parab olic version of the results in CIJ. The same correctors 
constructed in Cu\ play an important role in the parabolic equation. On 
the other hand, when we consider the following Porous M ediu m Equation 
in perforated domain, the viscosity method considered in CLl ] can not be 
applied directly. The equation will be formulated in the following form: find 
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the viscosity solution u^{x,t) s.t. 

An- - dtu, = in Qt,< (= ^at x (0, T]) 

Ue = on 9zQT,a* (= on,* X (0, T]) (PMEl) 

Ue = ge on fia* X {0} 

with 1 < m < oo and a compatible which will be defined at Section 

m The Dirichlet boundary condition can be considered an obstacle problem 
where the obstacle imposes the value of the solution is zero in the peri- 
odic holes. And the diffusion coefficient of [PME^) is mu^~^ and will be 
zero on dVLa* , which makes important ingredients of the viscosity method for 
uniformly elliptic and parabolic equations inapplicable without serious mod- 
ification. Such ingredients will be correctors, Harnack inequality, discrete 
gradient estimate, and the concept of convergence. Therefore the control of 
the degeneracy of {PMEl) is a crucial part of this paper. 

One of the important observation is that Ue{x,t) = '^^ will be a 

(i+t)™-i 

self-similar solution of [PMEl) if ^e{x) satisfies the nonlinear eigen value 
problem: 

A(/9e + =0 lu Via* 

= on dVLa* . 

The equation for ip^ is u niformly elliptic with nonlinear reaction term. The 
viscosity method in CL| . can be applied to the homogenization of ^p^ with 
some modification because of the nonlinearity of the reaction term. It is cru- 
cial to capture the geometric shape of p^ saying that ife is almost Lipschitz 
function with spikes similar to the fundamental solution of the Laplace equa- 
tion in a very small neighborhood of the holes. It is not clear whether we can 
find the geometric shape of ^p^ if we construct p^ by the energy method since 
i^^-weak solutions may have poor shapes. And then such self-similar solu- 
tion, U^{x,t) will be used to construct super- and sub-solution of (PMEl) in 
order to control the solution, u^{x,t), especially the decay rate of u^{x,t) as 
X approaches to the holes, T^^ at Section |H With the help of such control, we 
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are able to prove the discrete gradient estimate of the in order to compare 
the values of on a discrete lattice created periodically by a point in a cell. 
And we also able to show the almost flatness saying that the values of at 
any two points in each small cell are close to each other with an e-error if 
those points are away from the very small neighborhood of the hole in the 
cell. 

It is noticeable that the homogenized equation is expressed as a sum 
between the original equation and a term depending on the capacity and {if — 
u)+ as the case in the heat equation, Theorem 12.31 We also prove that such 
decoupling of terms will happen in the homogenization of porous medium 
equations in perforated domain, Theorem 13.61 But it is not clear whether 
such decoupling property holds in the general fully nonlinear equations of 
non-divergence type. Theorem 14.81 

This pa per is divided into three part: In Section m we review some fact 

ating obstacle problem for Harmonic operator) 



studied in 



CL| (highly oscil 



and extend the results of |CLj to the Heat operator and fully nonlinear 
parabolic operator. In Section [3l we study the elliptic eigenvalue problem 
in perforated domains, which describe the behavior of solution of porous 
medium equations at a point close to the boundary. And, in Section HI we 
deal with the estimates for the porous medium equation in fixed perforated 
domain. 

Notations: Before we explain the main ideas of the paper, let us summarize 
the notations and definitions that we will be used. 

• Qt = n X {0,T], diQr = dn x (0, T] 



• Ta^, M^^, fla^ and QT,a^ are described in Section [TJ 



• We denote by the cube {x = (xi, • • ■ , Xn) G M" : |xj — mj| < | (z = 
1, ■ • ■ , n)} where m = (mi, ■ ■ ■ , m„) G M". 
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Denoting by the corrector described in Section 2 in CL |. 



Numbers: a* 



for n > 3 and a* = e ^ for n = 2. 

n—2 — e 



2. Highly Oscillating Obstacle Problems 



Firs t, w e review results on the correctors in CL| . Likewise Laplace oper 



ator in 



CL| . the correctors will be used to correct a limit u{x, t) of a solution 



u^{x,t) for the obstacle problems i \H^\) in this section. 

Any possible limit, u{x,t), can be corrected to be a solution of each e- 
problem, ( |ij"eD , and it is also expected to satisfy a homogenized equation. 
The homogenized equation comes from a condition under which u can be 
corrected to u^. To have a oscillating corrector, let us consider a family of 
functions, We{x), which satisfy 




m 

Lie 

in Ta, 



(2.1) 



for some k > 0. In CLj, Caffarelli and one of authors construct the super- 
or sub-solutions through which they find the limit of depending on the 
decay rate of the size of the oscillating obstacles, a^. 

The next Lemma tells us that there is a critical rate a* so that we get a 
nontrivial limit of correctors, We- The reader can easily check following the 
details in the proof of the Lemma 2.1 in CLj. 



Lemma 2.1. Let = Coe". There is a unique number a* = 



lim inf 
lim inf 
lim inf 



-oo, 




1 



for any k > if a > a^, 
for a = and k = cap{Bi] 
for any k > if a < a*. 



In addition, we can also obtain the interesting property from JKO[ |. 
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Lemma 2.2. Set a = a^, then the junction satisfying 



Wf = I — Wf 



converges weakly to 1 in Hl^^{W^). 

2.1. Heat Operator 

We are interested in the limit u of the viscosity solution of {H^) as e 
goes to zero and the homogenized equation satisfied by the limit u. As we 
discussed in the introduction, there will three possible cases depending on 
the decay rate of a^. 

Theorem 2.3. 



Let u^{x,t) be the least viscosity super solution of\H, 



1. There is a continuous function u such that u in Qt with respect 
to U'-norm, for p > 0. And for any 5 > 0, there is a subset Ds C Qt 
and €0 such that , for < e < €0, ^ u uniformly in Ds as e — )■ 
and \Qt\Ds\ < S. 

2. Let a* = e"* for a* = ^^^2 /^'^ n > 3 and a* = for n = 2. 

(a) For Cott* < < Coa*, u is a viscosity solution of 

H[u]+kB^^{(p - u)+ = in Qt 

u = on OiQt. 

m(x,0) = g{x) onQ X {0} 

where ks^^ is the harmonic capacity of B,.^ if = lim£_^o ^exists. 

(b) If = o{a*) then u is a viscosity solution of 

H[u] =0 mQT 
u = on OiQt. 

u{x, 0) = g{x) on Q X {0} 
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(c) If a* = o(ae) then u is a least viscosity super solution of 



H[u] < 
u > (p 



in Qt 
in Qt 
on OiQt- 
on f2 X {0} 



u = 



u{x,0) = g{x) 



Remark 2.4. 1. The boundary data above can be replaced by any smooth 
function. And H[u] = f{x,t) can replaced by the heat equation. 

2. Ta^ = {a^x : x G D} can be any domain with continuous boundary 
as long as there is two balls B,.^ G D G Br^ for < ri < ri < oo. 
-Bri G D G Br2 is enough to construct super- and sub-solutions and 
then to find the behavior of correctors, Lemma \2. 1\ Then k = cap{D) 
and ks^^ = ko- 

2.2. Estimates and Convergence 

Every e-periodic function is constant on e-periodic lattice eZ". The first 
observation is that the difference quotient of u^, instead of the first derivative 
of tie, is uniformly bounded. The next important observation is that a suitable 
scaled is very close to a constant multiple of a fundamental solution in a 
neighborhood of the support of the oscillating obstacle, T^^ and that will 
be almost constant outside of it. These observations will be proved in the 
following lemmas. 

2.2.1. Estimates of 

Lemma 2.5. For each unit direction e G Z", set 



AX(x,t) 



Me(x + ee, t) — t) 



e 



Then 



/S,lu,{x,t)\ < C 



uniformly. 
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Proof, can be approximated by the solutions, m^^^, of the following penal- 



ized equations, [Fr|, 



-H[u,^s] {x, t) + I3s{u,^5{x, t) - (p^{x, t)) = in 

M,,5(x,t) = ondiQr (2.2) 

Me,5(x,0) = g{x) on X {0} 

where the penalty term /^^(s) satisfies 

;5Ks)>o, /3;(s)<o, ;g5(o) = -i, 

/^^(s) =0 for s > 5, I3s{s) -oo for s < 0. 

Let Z = sup(^ j)gQy |AgMe,5p and assume that the maximum Z is achieved at 
(xo,to)- Then we have , at {xo,to), 

H[\Alu,^s\^] < 0, and V|AX,5|' = 0. 

By taking a difference quotient, we have 

-i/[AX5] +/3K-)(AX^(x,t) - Al^.ix,t)) = 0. 

Hence 



-i/[|AX5|2] + 2|V(A>,/'2 



5. . 

2 



+ 2/3K-)(l AX5(a;, W - Kueix, t)Al^,{x, t)) = 0. 

Since the set Ta^ is e— periodic and is C^, we know |Ag(/?e| < C uniformly. 

If Z = |Agtie^5p > |Ag(y9ep at an interior point {xQ,to), we can get a 
contradiction. Therefore Z < lA^ip^l"^ in the interior of Qt- On the other 
hand, > ip and then Ps{ue,s — fe) = on a uniform neighborhood of 
OiQt- From the C^-estimate of the solution for the heat equation, we have 
|AgMe,5p < C supg^ < C* supg^ \ip\ on OiQt- Hence, by the maximum 
principle, Z < C{\\ip\\c^Qj,) + \\g\W{n))- □ 

Corollary 2.6. we have 
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\u^{xi, t) — Me(x2, t) \ < C{\xi — foT a UTiiform constant C when Xi — X2 G 

Lemma 2.7 (Regularity in Time). 

\\DMxM<C. 

Proof. Let T = sup(^ j)£Q^ |(^e,<5)tP and assume that the maximum T is 
achieved at Then we have , at 

H[\{u,^s)t\^]<^ 

By taking a time derivative in f l2.2p with respect to time t, we have 

+ /3K-)(K5)t - (^e)t) = 0. 

H6I1C6 

-i/[|(M*l']+2|V(Mtl' 

+ 2/3K-)(l(M*l'-K5)*(^.)*) = o. 

We also know |(</?e)i| < C uniformly. 

If T = > at an interior point (xi,ti), we can get a 

contradiction. Therefore T < < C for some C > in the interior of 

Qt- On the other hand, = u^^s > ip^ on diQr and ^^^^(x, 0) = (^(x) > ip{x, 0) 
then /^^(Me,^— v^e) = on a small neighborhood of dpQx- By the C^-estimate of 
the solution of the heat equation, we get the desired bound on the boundary 
and the lemma follows. □ 

Lemma 2.8. When a < a*, satisfies 

ip{x, t) - C—^ < t) 

for some /3 > n. In addition, there is a Lipschitz function u, such that 
1. 

~C g_2 < Ue{x,t) — U < 

, which implies the uniform convergence of to u. 
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2. u is a least super- solution of l[2^) in Theorem Ii2.3\) . 



Proof. (l.)Since > f m Ta^ x (0,T], we show that the inequahty can be 
satisfied in (5r,a,- For a given 6q > 0, let 



h^{x) = k sup 



Then H[h^] > c^k for (5 > n and a uniform constant cq. In addition, > 
he > in ^e- For any point (xq, to) in Qx^ae-, we choose a number po and 

large numbers /c, M > such that 

h{x,t) = - a^oP + Vv2(xo,to) ■ (a^ - Xq) 

+ ip{xo, to) + K{x) + M{t - to) < < t) 
on dBp^,{xo) X [ito,to] and Bp„(xo) x {^to} and 

h{x,t) < (p{x,t) < Ue{x,t) 
on {dTa^ n i?pQ(xo)} X [^tojto]- By the choice of numbers, we get 

H[h] = ^-M> 0. 

Therefore h < in {i?pp(xo)\Ta^} x [^toito], which gives us ip{xo,to) — 
C-^ < Ue{xQ,to)- Moreover, the least super solution v{x,t) above ^{x,t) 
is greater than which is the least super-solution for a smaller obstacle. 
Similarly, the lower bound of above implies v{x,t) < + C-^. 

Since is a super-solution, this implies (2). 

□ 

Lemma 2.9. 

Seta, = (^)i/2. Then 

osc Ue = 0(e'^), 
for m G eZ" fl supp (p and for some < 7 < 1 . 
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Proof. If we make a scale v^{x, t) = u^{a^x + m, a^t + tg), we have a bounded 
caloric function in a large domain {B^i2aX^)\Ba*jaX^)^ ^ [O5 ^^/'^e] such that 
v^{x,t) > (p^{aeX + m,t) > ip{m,t) — Ca^ on -Ba*/a,(0) x [0,e^/a^]. We may 
expect almost Louville theorem saying that the oscillation on the uniformly 
bounded set is of order o{e'^). Let be a caloric replacement of in B^j2a^ x 
[0,e^/a ^]. Then oscB^x[e2/a2-i,e2/a2] w;^ = o(e''') by applying the oscillation 



lemma, 



LB| |. of the caloric functions inductively: 



OSC We < 5o OSC We 

SH(a:o)x[to--R2,tQ] S4^(a:o)x[to-(4fl)2,to] 

for some Q < 6q < 1 and 7 ~ log^ (y) = . It is noticeable that 

5o = 1 — ^ for Ci > which comes from the Harnack inequality, 

sup w <Ci inf w 

5fl2 , „2l Sjj(a:o)x[to-fl2,to] 



S 



j/2{^o)xIto-T->'o--R^l 



for a positive caloric function w. Then the error = — is also a caloric 
in {5e/2a.(0)\5,*/,^(0)} X [0,eVa2] and = on [dB^/^aX^)] x [0, eVa^] 
and i?e\ae(0) X {0}, we also have 

< f < 2 sup 

Qt 

in {-Be/2a, (0)\-Ba*/a, (0)} X [0,e^/a^]. Since the harmonic function can be 
considered a stationary caloric function, we have 

< v{x,t) < (supv?) *^""^ , r = \x\ 

which means oscaB^x[o,i] v = 0(e"~^). Therefore we know 

OSC Ue = OSC We = o(e'^). 

SSag (m) X [to-<i2,to] SSi (0) X [(e/ae )2 - 1 , (^/ae )2] 



□ 



By the Lemma 12.51 12.71 and 12. 9[ we get the following corollary. 
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Corollary 2.10. we have 

\u^{x,ti) - u,{y,t2)\ < Ci\x -y\ + Cal^i - ^2!^ + o{e^] 

for {xM)AyM)e (u„e,z9E,,Hnfi) x(o,t]). 

Lemma 2.11. 

Seta, = Then 

osc = o(?7(e)) 

{B, (■m)\Ba, (m)}x [to -a^ M 

for m G eZ" fl supp (f and for some function rj^e) satisfying 

77(e) as e 0. 



Proof. The Lemma l2^ tells us that is almost constant on a set {dBa^ (m)} x 
[to — CL^,to] whose radius is greater than a critical rate a*. Let 

Ue{x, t) = sup Uf. 
aBa,x{t} 

for e {Q^ n fi} X (0,T]. Then, by the corollary ElDl we have 
\u,{x,ti) -Ue{y,t2)\ < Ci\x - y\ + C2\ti - t2\^ + o(e^) 

and 

\u,{z,t) -u,{z,t)\ < Ce^ 

for all (x,ti),(y,t2) e X (0,T], G {U„g,z"<95,^ M H fi} x (0,T] and 

for some C < 00. Therefore there is a limit u{x,t) of u,{x,t) such that 

sup |Me(x, t) — ^(x, t)| = o(?7(e)). 

{U„g5zn (9Sa, (m)nC} X [0,00) 

for some function 77(e) which goes to zero as e — 0. This estimate says that 
the values of [u{x, t) — Cr]{e))xTa P^^J^ obstacle below with a slow 

decay rate, a, » a*, in Lemma [2.81 which will gives us the conclusion. □ 
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2.3. Homogenized Equations 

In this section, we are going to find homogenized equation satisfied by 
the hmit u of through viscosity methods. 

1 

Lemma 2.12. Let a* = e"* for a^, = forn > 3 and a* = e ~^ for n = 2. 

Then for CqQ* < < CqO*, u is a viscosity solution of 
( 

An + kb^^ (v? - u)+ - = m Qt 

< M = ondiQT (2.3) 

u = g{x) m f2 X {t = 0} 

where KBrQ is the capacity of Br^ if tq = hme_i.o ^ exists. 

Proof. First, we are going to show that n is a sub-solution. If not, there is a 
quadratic polynomial 

P{x,t) = -d{t - t°) + l-aij{xi - - x°) + bi{xi - x°) + c 

2 ■' 

touching u from above at {x^ ,t^) and 

H[P] + K{ip - P)+ < -^0 < 0. 

In a small neighborhood of By{x^) x [t^ — ?7^,t°], there is another 

quadratic polynomial Q{x,t) such that 

H[P] < H[Q] in X - t^] 

< g(xO,tO) < P(xO,tO) -5o 
Q(x, t) > P(x, t) on 9P^(x°) X [t" - 7]'^, t^] and x {t° - r/^}. 

In addition, we can choose an appropriate number eo > so that Q satisfies 

H[Q] + - + eo) := H[Q] + <o < -y < 

and 

\Qix, t) - Q{x\ t°) I + \ip{x, t) - I < eo 
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in _B^(x°,t°) X [t° — ?7^,t°]. Let us consider 

Qe{x, t) = Q{x, t) + We{x)^o- 

Then we have 

H[Q,{x,t)]<-i^<0 

and 

g,(x, t) = g(x, t) + t°) - t°) + eo) 

> Q(x, t) + t°) - Q{x\ t°) + eo) 

on {Ta^ nB,,(x°)} X [to — ?7^, to]. Hence, by the maximum principle, Qe{x, t) > 
iPeix,t) in X [t° - r/^t°]. 

Now we define the function 

mm(u^,Qf:) X e By{x'^) 
Ue X e 

Since (^^^ ^ = o(e) as e — )■ 0, by Lemma 12.111 converges uniformly to u 
in Q. Hence, for sufficiently small e > 0, > on dBn{x^) x [to — ''7^, to]. 
Thus the function is well-defined and will be a viscosity super-solution of 
(12.31) . Since Me is the smallest viscosity super-solution of (12. 3p . 

Ue < Ve < Qe- 

Letting e 0, we have M(x°,t°) < (5(x°,t°) < P(x°,t°) = M(x°,t°) which is 
a contradiction. By an argument similar to the proof of Lemma 4.1 in Cu\ . 
we can show that u is also a viscosity super-solution of (12.31) . □ 

Lemma 2.13. 

When = 0(e") for a a > a^, the limit u is a viscosity solution of 

^H[u]=0 mQr 
u{x, t) = on diQr 

u{x, 0) = g on Q X {0} 



17 



Proof. For e > 0, H[Uf:] < 0. Hence the limit also satisfies H[u] < in a 
viscosity sense. In order to show m is a sub-solution in Qt, let us assume that 
there is a point (xq, to) £ Qt such that H[P]{xo, to) < — <^o < for a quadratic 
polynomial P such that (P — u) has a minimum value zero at (xo,to)- We 
are going to choose a small neighborhood of {xo,to), -B^(xo) x [to — ?7^,to], 
and a quadratic polynomial Q{x,t) such that 



(5(x,t) > P{x,t) on dBr^{xo) x [to - r^^to] and 5r,(xo) x {to - V^} 
H[Q]>H[P] mB,{xo)x[to-v\to] 
Qixo,to) < P(xo,to) - So 



Let Qe = Q{x,t) + {we — mint^e). Then H[Q^] = and Qe > (fe since 
1 — min We — J- oo as e — !■ 0. Hence min{u^, Q^) is a super-solution of ( |ffeD ; but 
min(Me, (5£)(xo, to) < Ue(xo,to), which is a contradiction against the choice of 
u,. □ 



Lemma 2.14. 

When tte = 0(e") for a < a^, the limit u is a least viscosity super solution of 



H[u] < 
u = 

u > ip 

^u{x,0) = g{x) 



in Qt 
on OiQt 
in D 



Proof. The proof is very similar to that of the lemma 4.3 in CL| . Likewise 
we only need to show u > ip. Let us assume there is a point (xo,to) such 
that u(xo,to) < (p{xo,tQ). We are going to construct a corrector with an 
oscillation of order 1, which is impossible in case that the decay rate of is 
slow, Lemma [2.11 and 12^81 For small e > 0, we have 

|M,(a;,to) - u{x,to)\ < i|n(xo,to) - v5(a;o,to)| 
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on (i?^(xo) n fiaj X {to}- For a sufficiently large constant Mi, we set 

Me + Mi\x — Xol > ^{x, to) 
on dBn^Xo) x {to}- Then we can set a periodic function 



= mm 



■m) 



for a sufficiently large constant M2 > and then it is a super-solution such 
that max We— min We > j\u{xo, to)—(p{xo, to) | > for small e > on Br^{xo) x 
{to} . Hence we can extend periodically W^ so that we have global periodic 
super-solution. But will not go to as e — 0, which is a contradiction 
against Lemma (12. ip . □ 



Proof of Theorem [HI (1.) Set D = (Ue<e„ Um&eZ" tX ' ) n For any 
d > 0, there is e^, > such that \D\ < 6. Corollary 12.61 shows the uniform 
convergence of on Q\D. 

(2-a),(2-b), and (2-c) come from Lemma 12.121 , Lemma 12.141 and Lemma 
EH □ 

3. Elliptic Eigenvalue Problem in Perforated Domain 

Before we deal with e-problem for the porous medium equation, we con- 
sider nonlinear eigen value problem, which will describe the behavior of the 
solution for the porous medium equation in a neighborhood of dQ. Let's 
consider the solution fe{x) of 

Aip^ + = 0, < p < 1 in 

V^e > in (EVe) 

= on Ta* U dQat - 

3.1. Discrete Nondegeneracy 

We need to construct appropriate barrier functions to estimate the dis- 
crete gradient of a solution ip^ of dEV^P on the boundary. 
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Lemma 3.1. 

For each unit direction Ci and x G dVL, set 

K^^ = 1 

and 



\\Al^.{x)\\ = 

Then there exist suitable constants c > and C < oo such that 

c < \\Al^,{x)\\ < C 

uniformly. 

Proof. Let /i+ be a solution of 

A/i+ = -MP - 1 in 
h+ = on dn 

with M > supq, (p^. Then, we have 

ip^ < in 
by the maximum principle and 

sup |V/i+| < C 

by the standard elliptic regularity theory. Thus, for x G dfl, 

\\Al^,{x)\\<\\Vh+{x)\\<C 

when we extend </?e to zero in M"\f2. 

To get a lower bound, we first show that the limit function (/?, of Lp^, is 
not identically zero. Let 

K = min ||V(^e||j;^2. 

^,eHi(o,*),||^,|| =1 
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f2„» 



For < ?7 G C^{Q) and corrector given in Section m set 6{x) = 'r]{x){l — 
w^{x)). Then 

/ \Vri{l -w,)\'^dx = V[r/(1 -We)] ■ V[r/(1 -wjjrfx 

(1 - w^Y\Vr]\'^ - 2?7(1 - w,)Vr] ■ Vw, + T]'^\Vw,\'^dx 
<2 {1 -w,y\\/r]\'^ + r]'^\Vw,\'^dx. 
Since ^ iVw^pfia; < C for some < C < oo, we get 

j \\/7]{l -W,)\^d^ <Ci 

for some constant < Ci < oo. On the other hand, 

for some constant C2,e depending on e. Since {1 — w^) ^ \ m. L'^iVL). we get 

hmCof = C2 < +00. 

Thus ^ 

1 \ ~^ 

7]{l-W,)\P^^dx 



and 

2 



Therefore we have 

Ae < — < +00. 

Then, the sequence {e} has a subsequence which we still denote by {e} such 
that 

°^ ^ (3.1) 

A, ^ A. 
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Since Hq is compactly embedded in L^"*"^, ip^ ^ ipin Lp~^^{Q) implies ||v5||lp+i(q) 
1. Note that A 7^ 0. Otherwise (f satisfies 

A<p = inn 
(f =0 on dQ. 

Then ip = which gives a contradiction since ||<^||Lp+i(n) = 1- For each A^, 
1 

the function ip^ = Xl~^ (p^ can be the solution of ( |EVJ . By fl3.ip . 

ip^ = Xl'" ip^ X^ip = ip m Hl{Vl). 

Since ||</3||lp+i(q) = A~ ||<^||lp+i(q) > 0, there is some constant > such 
that 

> 5o > mD dVt and ^ 0. 

Now we consider the ^eXo ~ denote by the minimizer of 



in K, = {'0, G H^i^e), i^e > ^e}- Thcu, by the Theorem 3.21 in |,TKOj . 
ip^ ^ ip in Hq{Q) and 

iP>Tp = ipx^ > (5o in D 

AV^ - = in 9\D 

for K = cap(i?i). Since ip^ satisfies the harmonic equation in VLa*\{il)e = ^e} 
with ip^ = ip^ = {] on dQa*, we get (p^ > ip^ in Qa*\{'^e = then in a 

neighborhood of dQ by the maximum principle. On the other hand, by the 
Hopf principle, 

inf \Vtp{x)\ >6i>0. 

Hence, there is a lower bound of ||Ag'?/'e||(x) if a; G dQ, which means 

IIA^^J > ||A^^/^J > 5i > foTxedQ. 

Therefore, || A|(^e|| is bounded below by some constant and the lemma follows. 

□ 



22 



3.2. Discrete Gradient Estimate 

Lemma 3.2 (Discrete Gradient Estimate). For the solution ip^ of (EV^), 
for all X E Q when we extend ipe{x) to zero in M.'^\Q. 

Proof. Let Gq and Gn^a* are the Green functions of the Laplace equation in 
fl and fla*, respectively. We choose constant 7 such that 

B2j{y) c ^a* {y e f^^i) 

and let the function Gn,a*,7 to be a solution of 

AGn,a-,7 = infi,.\5^(y) 

Gn,a:,'r{x, y) = Gn{x, y) on B^{y). 

Then we get 

Go-^aia ^ in Vtai\B^{y). 
Therefore, we obtain 

\^lGn,aiA<\KGn\ on an. (3.2) 

To get the estimate on dB.y{y), consider the difference 

G(x, y) = Gn{x, y) - Gn,a*,7(a;, y). 

Then G{x,y) satisfies 

AG = in Qat\B,iy) 

G = ondnudB^iy) 
G{x,y) = Gn{x,y) on dTa*. 

Note that G^{x,y) has similar behaivour to 0{\x — ?/P~") as \x — y\ — 0. 
Thus, 

max Gn{x,y) < - min Gnix,y) 
23 



for a sufficiently small 7 > 0. Thus, there exists a constant C > such that 
G{x,y) < max Gn{x,y) < min Gnix,y) - max Gnix,y) 

< GTJx -y)=G - . ) on dT* . 

Thus 

G{x,y) < CV^{x-y) in Vlai\B^{y). 
Since r^(x — y) = Gn{x, y) = on dB.y{y), we have 

\^lGn,ata\ < \^lGn\ + |A:G| < \AlGn\ + lA^r^l on dB^{y). (3.3) 

To show the estimate at the interior points, we use the approximation method. 
As in |Fr|, Gr2,a*,7 can be approximated by the solutions, Gn,a| ,7,155 of the fol- 
lowing penalized equations, 

^Gn,ai,^,s + - Gn,at,y,s + Gn ■ C{x)) = in Q\B^iy) 

Gn,at,'y,5 = ondQ (3.4) 

GQ,a*,7,5(a;, y) = Gn{x, y) on B^{y) 

where I3s{s) satisfies 

/3Ks)>0, /35(s)<0, /35(0) = -l, 
/^^(s) = for s >5, hm/35(s) — ;> —00 for s < 

5-s-O 

and a e-periodic function ^(x) G C°° satisfies 

eri--2 



0. 



Similar to the proof of Lemma 12.2. ![ we get 



/\{\AlGn,atn,s?) -2\V{AlGnKnM 



(3.5) 
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Suppose that ma,XQ\B^(^y) | AgGo^aj^^^^l^ occurs at an interior point Xq. If 
\AlGn,ai,'y,sW^o) > \AlGnWxo), then by we get 

> A{\AlGn,at,,A') - 2|V(A:Go,a.,7,5)r 



Therefore 



\AlGn,ai,^,s\^ < lA^G^r in the interior of n\B^{y). (3.6) 



By dMl), O and ([MD 

\AlGn,at,^,5\ < \AlGn\ + \AIT^\ in n\B^{y). 
By taking 6-^0 and 7 — )■ 0, we obtain 



Since 



we get 



\AlGn,at\ < \AlGn\ + \Air^\ in 



AlGn^at{x,y) ■ ^l{y)dy 



< / {\AlGn{x,y)\ + \AlT,\)^l{y)dy. 



Since hm.^o l^e^-^f^l = iVgCnl and \VGn{x,y)\ ^ 0(|x-y|^ ") as \x-y\ 
0, we get, sufficiently small e > 0, 



|A>.(x)| < C / {\x-y\'-^ + l)ip^Mdy < 00 



and lemma follows. 



□ 
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3.3. Almost Flatness 

Lemma 3.3. Set a, = {'-^f/^. Then 

osc ip^ = o(e'^) 

dBa,,{m) 

for m G eZ" fl Q and for some < 7 < 1 . 

Proof. If we consider the scaled function Ve{x) = (p^{a^x + m), will be 
bounded in B^i2a\Ba*jat and = on dBa*/a^. also satisfies 

Aw, = -alvl 

Let be a harm onic replacement of in B^/2a,\Ba*/ac- Following the proof 
of Lemma 3.5 in CL| which is similar to Lemma we get osc^Bi Qe = o(e^). 
Let's consider 

h(r) = — - - 

with r = \x\ and M = supw,. Then 

Ah = -alM^ < -ay, = A{v, - g,) in B,,2a\Ba*ja. 

h>0 = V,-g, on d{ B,/2a, \Ba* /a, } ■ 

By the maximum principle, we get 

ge<Ve<ge + h<g, + C 

for some C > on dBi. Thus 

oscu. < osc q. + Ce^ < oie^). 
dBi dBi - \ J 

If we rescale back to ip,, we can get the desired conclusion. □ 
Lemma 3.4. Set a, = (^)^/^ Then 

osc (fi, = o(e^) 

Be{m)\BaAm) 

for m G eZ" fl Q and for some < 7 < 1 . 
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Proof. By the Section 2.1 in CL| . there is a periodic corrector having 
properties 

Aw, = k and \1 - w,\ < Ce^~^ in M^^ 

for > and n < (3 < 2{n — 1). Let L',N' > be the constants to be 
determined later. We define the barrier function 

w^(x) = [l - We{x)] + L'\x - + M' + N'e 

with 

M' = sup ip,. 

dBa, (m) 

Then we can select sufficiently large numbers k,L' ^ 1 and A^' so that We 
satisfies 

Awe < Aip, in 

We > (fie on dfla^. 

By the comparison principle, we get 

V?, < We in Qa, 

Similarly, we get 

We < (fie in fla^ 

where 

We = [we{x) — l] — l'\x — + m' — n'e, m' = inf (pe 

for sufficiently large numbers l',n' > 0. Since B^{m)\Ba^[m) is a small 
region, we get 

\we — We\ < o{e^) in B^{m)\BaXm) 

with < 7 < 1 and lemma follows. □ 
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3.4- Correctibility Condition I 

Likewise highly oscillating obstacle problems, we need an appropriate cor- 
rector. However, unlike the highly oscillating obstacle problem, the corrector 
should be a super-harmonic with 



w, = 



in T„ 



Since the solution ty^ of ( 12. ip is sub-harmonic with = 1 in T^* , it is natural 
to consider the function 

We = b — bwe 

for some constant b > 0. 

Lemma 3.5. 

Let kb e be such that 



A (6 - bw,) + {b- bw. 



Then we have 



—bn. 



-bAwe + ib- bw,y = kb,e- 



h-bP 



where = lim^^o kb e O'^d the harmonic capacity of 

Proof. Set Ve{x) = We{a*x + m), then satisfies 

-bAve + (a:)2(6 - bveY = hM? 



'e = ■ Vf e I = 

Thus, we have 



AVfdx = (a 



*\2 



in \B, 
on dBi 
on OQq" . 



kb,e - hP{l - Vefdx. 



Qn" \Bl 



28 



On the other hand, from the elhptic uniform estimates, GT|, v^^ v con- 
verges to a potential function v of Bi in C^-norm on any bounded set where 
An = on ]R"\i?i, v = 1 on dBi, and n — )■ as |a;| — t- oo. Then we get 

—b / AiLdx = —b Vff ■ vda^ 



= -b Vv^ ■ {-u)dax -b Vv ■ {-i')da, 

dBi -J dB^ 

as e goes to zero. And we also have 



- bK,^ = hm [ih,, - W){a:f ■ i^r] = -4z-Ah - ff) (3.7) 

where is the harmonic capacity of Bi and kb = lim£_>.o kh^e- If we multiply 
equation (13. 7p by Tq"^, we obtain 

—bK„ = kh — b^ 

where k^^^ is the harmonic capacity of Brg. □ 

3.5. Homogenized Equation 
Theorem 3.6. 

1. (The Concept of Convergence) 

There is a continuous function ip such that (p^ ^ ip in Q with respect 
to L'^-norm, for q > 0. and for any 5 > 0, there is a subset Ds G Q 
and Co such that, for < e < eo, p^ uniformly in Ds as e — and 
\n\Ds\ < 6. 



_ 1 

2. Let a* = e"* for a* = for n > 3 and a* = e ^ for n = 2. Then 
for CqQ* < < C^a*, u is a viscosity solution of 

A(/? — K^^y (y9 + (yj^ = m f2 

< (^ = ond^ (3.8) 

(y9 > m 
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where is the capacity of if = liiiie^o % exists. 



Proof. By an argument similar to the proof of Theorem I2.3[ (1) holds. 
(2). For e > 0, 

Hence, the limit also satisfies 

A<^-K^,,^V^ + V^^ < 

in a viscosity sense. Thus, we are going to show that (/) is a sub-solution. Let 
us assume that there is a parabola P touching u from above at Xq and 

AP - K„ P + PP < -26q < 0. 

In a small neighborhood of Xq, Bj^{xq), there is another parabola Q such that 

D'^Q > D^P in B^{xo) 

Q{xo) < P{xo) - 6 

Q{x) > P{x) on (95^ (xo). 

In addition, for Q{xi) = mmB,^(xo) Q{x), 

AQ - K^^^Qixi) + QixiY < -6o < and \Q{x) ~ Q{xi)\ < Cr] 
in Bfj^xo). Then the function 

Qe{x) = Q{x) - w^{x)Q{xi) 

satisfies 

AQ, + < AQ - Q{x,)Aw, + (1 - w^YQixiY + {C7]f 
in i?^(xo) n Vta^. By correctibility condition I, Lemma [3.5[ 



AQe + Ql<AQ + h 



Q{xi),e 



5n 



< AQ + kQ^,,) + ^ + {Cvr 
<AQ- ^,^^Q{x,) + Q{x,y + I < -| < 
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for small e,?7 > 0. So AQ^ + Qp < and > on (9{5p(xo) fl fi^J for 
some p > 0. By a comparison principle , we get 

Ue < Qe 

in -Bp(xo) n Qa^. Thus we get Qel^^o) > V^e(2;o) and then Q(a:o) > ^p{xq). 
On the other hand, Q{xq) < P{xq) — 6 < ^pIxq), which is a contradiction. 
Therefore is a viscosity solution of f l3.8p . □ 

4. Porous Medium Equations in a fixed Perforated Domain 

Now we can consider the following porous medium equations. The main 
question is to find the viscosity solution Ue{x,t) s.t. 



Am™ — dtUe 
u, = 

Ue = 9e 







in QT,at (= f^a* X (0,T]) 

on dlQT,at (= d^a^^ X (0,T]) 
on X {0} 



[PME]] 



where 1 < m < oo and g^ix) = g{x)C,{x) for a smooth function g{x) G 
satisfying 

< (5o < |v^| <c on an 

and a e-periodic function ^(x) G satisfying 



e = l inMVi, 



0<e<l, e = inT,., 
Ae = inM^\{T,* UM"^}, < 0, 

Set Vf = uT" which is a flux. Then f ^ satisfies 





= 0. 



(4.1) 



Ve " Af e — dtVe 
Ve = 

Ve = gT 



in Qt,. 
on SjQt,* 
on X {0}. 



{PMEl 



In this section, we deal with the properties and homogenization for the so- 
lution V. 
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4.1. Discrete Nondegeneracy 

Let ip^ be a solution of the boundary value problem 

j_ 

A(pe + ipr =0 in Via* 

= on dVLa* . 

It is easy to see that the function 

(A + t)"-! V"^- 
satisfies the equation 

\/'A"'^AK,A-(K,A)t = 0. 

As in Lemma [231 the rescaled function ip^^a^^x + m) approach the harmonic 
function in B_L.\Ba^ as e — j- 0. Hence, for sufficiently small e > 0, 
becomes almost harmonic near T^*. Thus, ip^ is equvalent to the e-periodic 
function ^ in (14. ID near T^. , i.e., there exist some constants < c < C < oo 
such that 

c^e ^ i ^ Cfpe near T^*. 
Therefore, we can take constants < A2 < Ai < oo such that 



for the solution of the initial value problem [PMEf). Therefore, by the 
nondegeneracy of yj^ in a neighborhood of dVt, we can get the following result. 

Lemma 4.1. 

For each unit direction e and x G dVt, set 

Ve{x + ee, t) — Ve{x, t) 



Then there exist suitable constants c > and C < 00 such that 

c < \ Alv^{x,t)\ < C 

uniformly. 



32 



4-2. Almost Flatness 

For small 60 > 0, we consider the set 

Ta* = {Ve < So}. 
As we mentioned above, satisfies 

K,A: < Ve < K,A2 

for some < A2 < Ai < 00. Hence the function is trapped in between 
amd V(;,A2 near the T^^ . Thus, there exists a uniform constant c > 1 such that 

Ta* C Ta* C Tea*- 

Therefore, the hole 7^* is not much different from Ta*. Since satisfies 

0<c<v,<C<oo, in flat X (0, T], (4.2) 

by uniformly ellipticity of v^, (14.21) has the harnack type inequality. Following 
the same argument in Section l2.1( Heat Operator), we have the following 
Lemma. 

Lemma 4.2. Set a, = (^)^/^ Then 

OSC = o(€'^) 

{S,(m)\B„,(m)}x[to-ai,to] 

for m G eZ n supp ip and for some < 7 < 1. 

4-3. Discrete Gradient Estimate 

can be approximated by the solutions, u^^^, of the following penalized 



equations, Frj, for sufficiently large number M > 0, 



Ai;,,5 - v'^^s ^{v,^5)t + P5{-v,,5 + 5 + M^{x)) = in Qr 
Ve,5 = S on diQr 
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where (3s{s) satisfies 

Ps{s)>0, /3^'(s)<0, f3s{0) = -l, 
l3s{s) = for s > (5, 
\im Bs(s) — — oo for s < 0. 

Using this, we obtain the following results. 

Lemma 4.3. // {v^)t is non-positive at t = 0, then {v^^s)t < for all t G 
(0,T]. 

Proof. First, we assume 

K,)t(-,0) <0. 

Since v^^s is positive, we have 
Hence 

Let fs{s) be a function having the maximum value of {v^^s)t at t = s, then 
there exist points x{s) = (xi(s), ■ ■ ■ ,x„(s)) G M" such that 

fs{s) = {ve,5)t{x{s),s). 

Since x{s) are maximum points, we have 

{f5)s = iiv,^s)t)s = iiv,^s)t)t + V(t;e,5)t ■ = {iv,^s)t)t- 

Hence 

(^)s < f 1 - ^ - vlrfsP's < C,,sfs 

, which implies 

< /5(0)e^-^^ < 0. 
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When (fe,<5)t(-, 0) < 0, we can approximate (fe,5)(-,0) by a smooth initial 
data, (fe,5,fc)(-, 0) such that (fe.^.A.Otl', 0) < 0. By the argument above, we 
know that {ve,5,k)t{-, s) < and then {v^,5,k){-, Si) > {ve,5,k){-, S2) for Si > S2. 
Since the operator is uniformly elliptic on each compact subset D of ^2^, we 
have uniform convergence of v^^s,k to v^^s- Hence Si) > {v^^s){-,S2) for 

Si > S2 and lemma follows. □ 

Lemma 4.4. // ^ | is non-positive, then 

\^Ve,s\L°° < 

with Ce satisfying Mm^^oCe = 00. 
Proof. For i G {1, ■ ■ ■ ,n}, we will have 



m 



Hence 

- Ks'\(^^,s)lX - 2/3'(-)(|(M.J' - KsUM^x.) > 0. 

Let Xi = sup(2,^^)gQ^ |(^e,(5)a;, P and assume that the maximum Xi is achieved 
at (xo,to). Then we have , at {xo,to), 

A(n,,5)' < and {{u,,5)i)t > 0. 



By the Lemma 14.31 we get 

Ks)t < 0. 

Thus, if Xi = Iv^^sll^ > l^x, P at an interior point (xcto), we can get a 
contradiction. Therefore Xi < in the interior of Qt- To get a bound 
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of the maximum Xi on the lateral boundary diQx or at the initial time, we 
consider the least super-solution / of the obstacle problem 

A/ < inVL 

f{x) > g{x) in n 

f{x) = 6 on d^l. 

Then / is a stationary super-solution with f > g in Q and / = Ve,s on dQ. 
Hence, by the maximum principle and Hopf principle, we get 

X^<C{U\\cHQr) + \\9\\cHQr) + \\f\\cHQT)) 

and the lemma follows. □ 

Lemma 4.5. For each unit direction e, we define the difference quotient of 
at X in the direction e by 

Ve,s{x + ee,t) - Ve,s{x,t) 



dt \t=0 



//"TrL-n non-positive, then 



\KVe,s\ < C 

uniformly in Qt- 

Proof. Since M^{x) is e-periodic, we will have 

A(A:K,)) - Al{v%-'){v,,s)t - v%~\KKs)), - KKs)P'{-) = 0. 
Hence 

a(|a:k,)P) -2a:k,)a:(4-')k,), 



V' 



,5 



(|A:K,)P)^-2|A:K,)P/3'(.)>0. 



——1 

Since Al(ve,s) and Ag(?;™^ ) have different sign, we can get a contradiction 
if |Ag(f£5)p has a maximum value in the interior. Hence, 

|A>,|2 < C int{QT) 

for some constant C > 0. On the lateral boundary, the estimate is obtained 
from the Lemma [4.11 Thus we get lAgW^I < C in Qt- □ 
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Corollary 4.6. If is non-positive, then we have 

{Ve)t < 

and 

|A>,| < C 

uniformly in Qt- 

Proof. By Lemma I4.4[ for each e > 0, v^^s converges uniformly to up to 
subsequence. Then v^^s{.x,ti) > v^s{,x,t2) for ti < t2 implies v^{x,ti) > 
v^{x,t2) and then {v^)t{x,t) < 0. By the Lemma [4.5[ we have 

v^^sjx + ee,t) - v^^s{x,t) ^ ^ 
e 

Therefore, by taking 6-^0, \Alv^\ < C. □ 

4.4- Correctibility Condition II 

Likewise elliptic eigenvalue problem, we need an appropriate corrector. 
Similar to the correctibility condition 1, we start with the following form 

We = d — dWe 

where is given by (12. ip and for some constant d > 0. 

Lemma 4.7. Let kc,d,e be such that 

((1 - w,y - l)dPc +{d- dw,YA{d - dw,) 

= ((1 - w,y - l)dPc - d'+^'il - w^YAw, = 

for some c,d > 0. Then, we have 

—d '^^Kg^^ = kc^d 

where k^^d = hm^^o kc,d,t one? is the harmonic capacity of B^q ■ 
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Proof. Set v^{x) = w^{a*x + m), then satisfies 

' ((1 - w,)P - l)dPc - d'+^il - v^YAv, = K^dM? in Qf\Br 
v^ = 1 on dBi 

=\i^ ■ Vfel =0 on OQq' . 

Thus we get 

{alYdFc j {{l-v,Y-l)dx = 

d'^P j (1 - v^YAvJx + (a:)2 J hAedx. 



Similar to the correctibihty condition I, Lemma [3.5[ letting e — )■ 0, we get 

- d'^'^n,^ = Ihn [kM<n^r] = (4.4) 

where is the harmonic capacity of Bi and kc,d = linie-j.o ^c,d,e since = 
(1 — We) ^ 1 in If we multiply equation (14. 4p by rg~^, we obtain 

where k^^ is the harmonic capacity of Bj.g. □ 

4-5. Homogenized Equation 

Finally, we show the homogenized equation satisfied by the limit u of 
through viscosity methods. 

_ 1 

Theorem 4.8. Let a* = e°* for = forn > 3 and a* = e ~^ forn = 2. 
Then for cqq* < < C^a*, v is a viscosity solution of 

;i~m(Af — Kg^^v^) — vt = in Qt 

f = on OiQt 

v = g^ zn 1] X {t = 0} 

where the capacity of B^^ if Tq = lim^^o ^ exists. 
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Proof. For e > 0, 



1--!- 



K„ Ve < 0. 



Thus, the hmit v of also satisfies 

v'^-^{Av - -vt<0 

in a viscosity sense. So we are going to show that i; is a sub-solution. Let us 
assume that there is a parabola P touching v from above at Xq and 

P^~^{AP - K^^^^P) -Pt< -26o < 0. 

In a small neighborhood of Xq, 5^(xo) x [tg — T]'^,to], we can choose another 
parabola Q such that 

in Bnixo) x [to - ?7^to] 
in Briixo) X [to - ?7^to] 



Qixo,to) < P{xo,to) - 6 



Q{x,t) > P{x,t) 



on 



{dB^{xo) X [to -^^^o]} 
n{5^(xo) X {to - 



and 



"(Ag-ztgi) -g^ < -(5o < 

for Qi = Q{xi,ti) = mmB^^xo)x[toWM Q{x,t)- Let us consider 

Qe{x, t) = Q{X, t) - QiWe{x) + €0 + h{x, t) 

for a small number < eo < | and a function h(x, t) we choose later. In 
{Br,{xo) n QaA X [to - ^7^to], Qe Satisfies 

ge"-Ag,-(g,)i< 



gi "^il-wX~-AQ-Q, '-{l-WeY-^Awe-Qt 



+ 



c{Q -Qi + eo + hy-^AQ + (g - Qiw, + eo + /i)^'^ A/i 



-h. 



[1] + [2] 
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with c = if AQ < and c = 1 if AQ > 0. To remove the [2], we consider 
the following initial value problem 

{a'^{x,t)Dijh -ht = f{x,t) in M" x (0, oo) 

h>0 in M" X (0, oo) 

h{x, 0) = Qiw^{x) 

with 

if ^ ^ J 

I- --,1 L ' 

[{Q - Qi)C{x, t) + eo + h\ otherwise 

fix, t) = -c[{Q - QiXix, t) + eo + h] AQ, 
C(x, t) e < C(x, t) < 1, Cix, t) = l in B,{xo) x [0, r]^] 

and C{x,t) = in{5,+^,2 X [0,(r/ + r/2)2]}c. 

Since the equation has non-degenerate coefficients, we can find the solution 
h{x,t) of the initial value problem. We can also observe the fact that the 
solution h{x, t) decays rapidly in a small time because ^ as e — j- in 
-ffQ(M"). Hence, for sufficiently small e > 0, we get 

^ h{x,t) < - at t = 7f. 

Therefore, satisfies 

qY^aq, - {Q,)t <q\~^aq + qY^ [(1 - w]~^^ - 1)] Ag 

-Q\~^{l~w,f'^^Aw,-Qt 
in {Bri^xo) n fiae} X [^0 — '7^5^o]- By correctibility condition II, Lemma W7l\ 

qI'- AQ, - {Q,)t <qI-AQ + - Qt 

<Ql-^AQ + k^^^^^ + ^-^-Qt 

< q'-- [AQ - ^,^^^Q^) + I - < -| < 
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1- — 

for small e > 0. Hence Qe " AQ^ — {Qt)t < and > on d{Bp{xQ) n 
^a,} X [^o~P^)^o] and {i?p(xo,to)nfia, } X {^o~P^} for some p > 0. By a com- 
parison principle, ^^(xcto) > ^e(2;o5^o) and then Q(xo,to) + f ^ ^(a^c^o)- 
On the other hand, Q(xo,to) < P(xo,to) — 5 < u(a;o,to) ^ <^0) which is a 
contradiction. □ 
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